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Abstract 

In this study we discuss the performance of approximate SQS supercell models in describing the cubic 
elastic properties of Bl (rocksalt) Tio.5Alo.5N alloy by using a symmetry based projection technique. We 
show on the example of Tio.5Alo.5N alloy, that this projection technique can be used to align the differently 
shaped and sized SQS structures for a comparison in modeling elasticity. Moreover, we focus to accurately 
determine the cubic elastic constants and Zener's type elastic anisotropy of Tio.5Alo.5N. Our best supercell 
model, that captures accurately both the randomness and cubic elastic symmetry, results in Cu = 447 GPa, 
C12 = 158 GPa and (^44 = 203 GPa with 3% of error and A = 1.40 for Zener's elastic anisotropy with 6% 
of error. In addition, we establish the general importance of selecting proper approximate SQS supercells 
with symmetry arguments to reliably model elasticity of alloys. In general, we suggest the calculation of 
nine elastic tensor elements - Cu, C22, ^33, C12, C13, C23, C44, C55 and Cqq, to evaluate and analyze the 
performance of SQS supercells in predicting elasticity of cubic alloys via projecting out the closest cubic 
approximate of the elastic tensor. The here described methodology is general enough to be applied in 
discussing elasticity of substitutional alloys with any symmetry and at arbitrary composition. 
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INTRODUCTION 



TiAlN coatings with their good oxidation resistance and excellent mechanical properties have at- 
tracted high technological and academic interest [Ij. Several studies have been devoted to discuss 
these alloys from different aspects to extend our understanding in maximizing their functionality and 
operational efficiency. The thermodynamics, phase stability and spinodal decomposition in TiAlN 
have been analyzed [21 |3j, also on the influence of nitrogen off-stoichiometry [3] and pressure [5]. 
Furthermore, the theoretical prediction of the mixing enthalpy when alloying TiAlN with Cr has 
resulted in a general design route to improve the thermal stability of hard coatings |B] • Recently, the 
importance of the significant elastic anisotropy in TiAlN on the isostructural spinodal decomposition 
has been disccussed [ZIIH]- Though, the available theoretical tools with the help of modern supercom- 
puters allows us to tackle such complex physical phenomena in alloys [9] , the prediction of anisotropic 
tensorial materials properties of substitutional alloys from first principles remains a challenging and 
highly requested task in computational materials science (TUl [H]. For example, in dynamical sim- 
ulations suitably designed simulation cells can greatly reduce the computational costs of predicting 
the temperature dependence of the elastic and piezoelectric tensors of alloys. The importance of the 
elastic and piezoelectric tensors of materials can be underlined not only by its fundamental role in 
materials science but also their distinguished usage in (micro) mechanical modeling, engineering or 
designing of machine elements, sensors, telecommunication devices, aircrafts, etc.. 

Although the ordinary scalar cluster expansion [12] offers an exact treatment of the thermody- 
namics of alloys and its tensorial generalization [TT] gives the most elegant description of anisotropic 
tensorial materials quantities of alloys, the computationally less demanding and less complex special 
quasirandom structure (SQS) approach p3] is more favorized due to its simplicity and success. For 
example the giant piezoelectric response of ScAlN alloys [HI [15] or the mechanical properties of TiAlN 
[7| have been successfully described within this approach. Using different superstructures, Mayrhofer 
et al. have discussed the impact of the microscopic configurational freedom on the structural, elastic 
properties and phase stability in TiAlN [16]. In B-doped wurtzite AIN significant configurational 
dependence of the piezoelectric constant has been predicted [17J with presuming wurtzite symmetry, 
similarly to the discussion of electronic properties and nonlinear macroscopic polarization in III-V 
nitride alloys [l8l[T9]. 

In fact, in these studies the success of the SQS approach in describing the energetics of alloys is 
presumed for predicting tensorial materials properties when using different approximate SQS super- 
cells or even ordered structures. These works were either only predictive on the materials constants 
or the confirmation of the applied approximate structural model was based on the experimental 
agreement of the results. Moreover, most of the previous theoretical works on predicting elasticity 
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and piezoelectricity of alloys presumed the experimentally observed symmetry for the modeling SQS 
supercells, though the substitutional disorder of the atoms in general breaks the local point symmetry 
of the supercell, and focused only on the corresponding principal symmetry non-equivalent tensor 
elements. While the symmetry arguments based tensorial version of the cluster expansion [H] gives 
an exact approach to completely include the local point symmetry of the materials, improperly cho- 
sen SQS supercells may result in large discrepancy between theory and experiments or in erroneous 
theoretical findings. 

The SQS approach, in principle, is not aimed to generate structures with the inclusion of local 
point symmetry and thus to provide the proper, full description of tensorial properties of alloys. In 
fact, different SQS supercells break the symmetry somewhat differently and thus the comparison of 
the differently shaped and sized SQS supercells in terms of modeling the elasticity of cubic Tio.5Alo.5N 
is a rather complex issue. Hence, detailed systematic studies on the application of SQS supercells in 
predicting elastic constants of alloys are required to establish their performance and to determine their 
applicability limits. For example, J. von Pezold et al. |2U] have recently evaluated the performance 
of symmetricaly shaped - {A x A x A), supercells for the description of elasticity in substitutional 
AlTi alloys and obtained convergence and error bars for the cubic-averaged principle cubic elastic 
constants within the supercell configuration space. However, a general concept of comparing and 
measuring different sized and shaped SQS supercells in describing tensorial materials properties is 
still lacking. 

In this study, we present a general projection approach to establish a way of comparing the 
ab-initio calculated elastic constants of Bl Tio.5Alo.5N obtained with applying different sized and 
shaped SQS supercells. We accurately predict and extensively discuss the calculation of the principal 
cubic elastic constants of Bl Tio.5Alo.5N within the SQS approach. In general, we establish the 
importance of selecting proper SQS supercells with symmetry arguments to reliably model elasticity 
of alloys. Namely, we show that supercells even with good short range order (SRO) parameters may 
result in large non-cubic elastic constants and, on the contrary, supercells with bad SRO parameters 
might approximate cubic elastic symmetry fairly accurately. We give the convergence of elasticity 
with respect to different SQS supercells for Bl Tio.5Alo.5N via the symmetry projected cubic elastic 
constants. Moreover, we suggest the calculation of 9 elastic tensor elements - Cu, C22, Cu, C13, 
C23, C*44, C55 and Cqq, instead of 21, to evaluate and analyze the performance of SQS supercells in 
predicting elasticity of cubic alloys. 
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METHOD 



In this section we provide a description of the techniques we apphed to calculate and analyze 
the approximate elastic constants of cubic Bl Tio.5Alo.5N alloy. First we explain the applied special 
quasirandom structure (SQS) approach of modeling alloys and discuss the difficulties of describing 
the proper symmetric tensorial materials constants within the model. After that, we summarize the 
computational details of obtaining the energetics and extracting the elastic tensors of the approximate 
SQS supercells. Finally, we present a general projection method that provides a technique to compare 
and analyze the calculated approximate elastic tensors and what establishes a principle to discuss the 
supercell models in terms of modeling elasticity with the inclusion of local point symmetry. The here 
described methodology is general enough to be applied for substitutional alloys with any symmetry 
and arbitrary composition. 

Special quasirandom structure approach and its symmetry 

The special quasirandom structure (SQS) approach [13] greatly reduces the computational difficul- 
ties of modeling thermodynamics, mechanical and electronic materials properties of random alloys. 
The approach models the substitutional disordered alloys with ordered superstructures. The basic 
structural element of the SQS model is a supercell, what is aimed to capture the structural short 
range order (SRO) in alloys while its periodic repetition introduces spatial long range order (LRO) 
[18]. The degree of SRO is usually measured by the Warren-Cowley parameter [21], which for a 
pseudobinary Ai-^^B^^N alloy is defined as aj = 1 — Pb{R)/xb, where Pb{R) is the probability of 
finding a B atom at a distance R from an A atom and xb stands for the concentration of B. A 
perfectly random alloy is characterized by vanishing SRO, while a > and a < define clustering 
and ordering, respectively. In terms of modeling disorder, approximate SQS supercells with with 
small or vanishing SROs up to a certain neighboring order can be compared if and only if the in- 
teraction parameters are also known. In this work the atomic configurations in the supercells were 
obtained by including the Warren-Cowley SRO parameters of the first seven nearest-neighboring 
shells. Namely, the disorder has been considered up to the seventh neighboring shell on the metal 
sublattice. Accordingly, the SRO parameters were calculated only on the Ti-sublattice. In order 
to achieve the closest possible model of the perfectly random alloy in the chosen sized and shaped 
supercell approximation {A x B x C) , a Metropolis- type simulated annealing algorithm [22] has been 
applied with a cost-function built from the properly weighted nearest-neighbor SRO parameters. 

The SQS supercell approach in general breaks the local point symmetry at different stages. The 
substitutional disorder changes the microscopic local environments which results also in some dis- 
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torsions on the lattice parameters. Namely, after full relaxation the supercells will have a general 
triclinic shape. Moreover, the SQS approach in modeling the substitutional disorder of alloys allows 
one to apply arbitrary supercell shape and size - {A x B x C) in terms of lattice vectors. This 
arbitrariness though increases the variational freedom to obtain closely vanishing SRO parameters 
with relative small supercell size at any alloy composition, it also spoils the symmetry of the model. 
Thus, the elasticity of the Bl Tio.5Alo.5N alloy is modeled with fully relaxed SQS supercells an it is 
described by 21 elastic constants, instead of the three principal cubic constants, Cn, C12 and C44. 
Namely, in the SQS approach the elastic tensor of the model belongs to a symmetry class that is 
lower than the one that the alloy shows experimentally. Furthermore, different SQS supercells break 
the symmetry somewhat differently, which means that the comparison of the results can only be 
done after certain alignment. In this study we show that a projection technique can provide such an 
alignment in the example of the Bl Tio.5Alo.5N alloy. 

Calculational technique to obtain the elastic tensors 

To obtain total energies and extract the elastic constants of the supercells introduced above. Den- 
sity Functional Theory (DFT) calculations have been performed with using the plane-wave ultrasoft 
pseudo-potential [23] based Quantum Espresso program package The exchange correlation en- 
ergy was approximated by the Perdew-Burke-Ernzerhof generalized gradient functional (PBE-GGA) 
[25] . The plane- wave cutoff energy together with the Monkhorst-Pack sampling [26] of the Brillouine 
zone were tested and sufficient convergence was achieved. The pseudopotentials were downloaded 
from the library linked to Quantum Espresso and tested by calculating the elasticity of bulk Bl AIN 
and TiN in agreement with literature values [TBI EZ] • In obtaining the ground state structure of the 
modeling supercells, both the lattice parameters and the internal atomic coordinates were relaxed by 
using the extended molecular dynamics method with variable cell shape introduced by Wentzcovitch 
[28] . Accordingly, during the relaxation the supercells geometries have been changed from the initial 
cubic-like lattice structure and converged to a slightly distorted triclinic shape with vanishing stress 
tensor. Thus, we avoid any residual structural stresses, which is essential in performing an accurate 
comparative analysis of the calculated elastic tensors. In this dynamics, a value of 0.02 KBar was 
taken as convergence threshold for the pressure. The elastic constants were calculated via the second 
order Taylor expansion coefficients of the total energy 

_ 1 g^E(6i,...,e6)) 
~ Vo de^de, 

where Voigt's notation is used to describe the strain e and elastic Cij tensor ji29l |30|. To obtain the 
entire elastic tensor, namely the 21 elastic constants of each supercell, 21 different distorsions have 
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been applied without volume conservation. The elastic constants were calculated by standard finite 
difference technique from total energy data obtained from ± 1% and ±2% distorsions. 

Projection of the elastic tensor to the closest elastic tensor of higher symmetry 

In this section we describe the projection technique introduced by Moakher et al. [52] to obtain 
the closest elastic tensor with higher symmetry class for any given elastic tensor with arbitrary 
symmetry. This projection technique allows us to extract the largest cubic part of the calculated 
elastic tensors. It introduces a tool to compare the obtained approximate elastic tensors and measure 
the appropriateness of the SQS supercells in modeling the elasticity of Bl Tio.5Alo.5N. 

The symmetric elastic tensor has 21 inequivalent elements for the most general triclinic system. 
A system with higher point symmetry requires less parameters in describing its elastic behavior. 
For example, with cubic symmetry the material has only 3 principal elastic constants, Cii,Ci2 and 
C44, while the hexagonal point symmetry results in 5 elastic constants, Cn, C12, C13, C33 and C44. 
Nevertheless, any elastic tensor can be expressed as a vector in a 21 dimensional vector space, with 
the following components 

E = (Cii, C22, C'sa, \/2C23, V^Cis, V^Ci2, 2C44, 2C55, 
2C66, 2Ci4, 2C25, 2C36, 2C34, 2Ci5, 2C26, 2C24, 

2C35, 2Ci6, 2a/2C56, 2a/2C46, 2A/2C45), (2) 

where the -\/2's ensure the invariance of the norm on the representation, whether it is vector or 
matrix. For the basis vectors see Ref.^T]. The following projectors Pgym generate the closest elastic 
tensor with higher symmetry via 

where -Esym has higher point symmetry. The term closest here is used in the sense, that the Euclidean 
distance \\E — E^'^'^W is minimum. 

To obtain the closest cubic approximate in our study of Bl Tio.5Alo.5N, we applied the projector 
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given as a 21x21 matrix, 
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Accordingly, the projected cubic elastic constants can be achieved via the following simple averaging, 

C12 

6 

C44 + C55 + Cee 

<-^44 = • [p) 

We can call them cubic-averaged elastic constants, since the equation is equivalent with averaging over 
the three orthogonal directions, [100], [010] and [001]. We note here that this averaging was used by 
von Pezold et al. in searching for optimized supercell in AlTi alloys. Thus, to obtain the closest cubic 
projection of an elastic tensor with arbitrary symmetry one needs to derive 9 different distorsion and 
calculate 9 independent tensor elements, like Cn, C22, 6*33, C23, C13, C12, C44, C55 and Cqq. In case of 
cubic symmetry Eq.(|5| results in the well-known cubic identities of the elastic constants, see Eq.([9]). 
For modeling elasticity in hexagonal alloys, one needs the closest hexagonal approximation that 
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can be obtained via the projector 
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that acts in the same 9 dimensional subspace and results in the following expressions for the projected 
hexagonal elastic constants, 

Cn = 3(Cn + C22)/8 + C12/4 + ^66/2, Cu = {Cn + C22)/2 + 3Ci2/4 - 

Ci3 = (C'13 + C23)/2, C33 = C33, C44 = (C44 + C^^)/2. (7) 

A detailed derivation of the projectors for the all symmetry classes, monoclinic, orthorombic, tetrag- 
onal, trigonal, hexagonal, cubic and isotropic can be found in Ref. [311 [32] . It is worth to mention 
that not all projectors can be defined in the above used 9 dimensional subspace. Furthermore, the ap- 
plication of this projection technique allows one to spilt the elastic tensor into a direct sum of tensors 
with different symmetry. Such decomposition is possible, for example, on the following routes, 

-^cubic -^'tetragonal ~t~ -^orthorombic 
~l~-£'monochnic ~l~ -^'trichnic 
-^hexagonal ~l~ -^tetragonal ~l~ -^orthorombic 
~l~ -^monoclinic ~l~ -^triclinic- 



(8) 



Thus, with calculating the norm of the components one gets information about the different contri- 
butions and can analyze elastic anisotropy in general |31j . 



RESULTS AND DISCUSSION 



In this section we present a comparative analysis of the calculated approximate elastic tensors 
obtained for the cubic (Bl) TiAo.5lo.5N alloy within the special quasirandom structure approach. To 
get different levels of the approximation of the elasticity in cubic Tio.5Alo.5N, several approximate 
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SQS supercell models have been generated with different shape and size, such as (2x2x2), (2x3x2), 
(4x3x2), (4x3x4), (4x4x3) and (4x4x4). Here SQS supercell sizes are measured in terms 
of the fee unit vectors. To have a more complete comparison of the calculated elastic tensors, we 
present results obtained with the ordered LIq structure and three other structures, denoted here by 
Cl-(2 X 2 X 2), C3-(2 x 2 x 2) and Bl-(2 x 2 x 2). These three structures are not based on the 
fee unit cell but on the fee Bravais cell. The Cl-(2 x 2 x 2) and C3-(2 x 2 x 2) was created by 
Mayrhofer et al. [16j with considering the number of bonds between the host and doping atoms. 
The C3-(2 x 2 x 2) structure was designed with preserving the cubic symmetry. The Bl-(2 x 2 x 2) 
structure was obtained by von Pezold with using a Monte-Carlo scheme and averaging over the 
three orthogonal main crystallographic directions. The SRO parameters of all superstructures are 
summarized in Table |lj In the case of the (2x2x2), (2x3x2) and (4x4x4) supercells, those 
atomic configrations have been chosen that resulted closest to randomness in our approximation, 
i.e. almost vanishing SRO parameters up to the seventh neighbor shell. The larger SROs in case 
of the (2x2x2) supercell are the consequence of the low configurational freedom in the supercell 
and indicate less perfection in the randomness. In the case of the (4x3x2) supercell size two 
different atomic configurations have been considered with very different SRO parameters. The * 
marks the SQS structure that is less random. The calculated SROs of the Cl-(2 x 2 x 2) and C3- 
(2x2x2) structures show alternating systematics that is related to the used construction strategy. 
The SRO parameters deviate considerably from zero in these two cases. For example, the cubic 
symmetric C3-(2 x 2 x 2) shows perfect ordering in every second neighboring shell. In comparing the 
SRO values in Table |l| the (4x4x4) supercell gives unambiguously the closest model of a totally 
random(pseudo-)binary alloy in our SQS approximation. 

The structural optimization of these supercells resulted in slight structural distorsions, what are 
summarized in Table [ITj Table [IT] gives the size resolved lattice parameters of Bl Tio.5Ao.5lN within 
the different SQS supercell models. The lattice parameters, especially the length of the cell edges, 
show some noticeable deviation from the cubic structure, but only for the LIq, (2x2x2) and Cl- 
(2x2x2) supercells. What correlates with the systematic alternation of the large SRO parameters 
of these cells. This suggests that the observed structural deviation is related to the low degree of 
freedom of internal atomic arrangement. In the other supercells with higher substitutional atomic 
disorder /randomness the cubic imperfection is nearly negligible. 

For each of these structures the full elastic tensor has been calculated. The elastic constants were 
obtained independently, as 21 different distorsions were applied. The obtained elastic tensors are 
summarized for all the structures in Appendix A. All the obtained tensors exhibit deviations from 
a strict cubic symmetry, in which the principal non-vanishing elements should show the following 
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relationships, 



tvll — t^ll — t^22 — t^33; tvl2 — tvl2 — <-^13 — <-'23, 

C44 = C44 = C55 = Cee- (9) 

Appendix E| also lists the elastic tensors of bulk Bl TiN and AIN obtained with the supercell size 
(4x4x3), where one finds the cubic symmetry of elasticit constants with numerical error. The values 
show good agreement with the literature data [27] obtained with different techniques. As the 
C3-(2 X 2 X 2) supercell preserves the cubic symmetry, its elastic tensor shows the cubic relationships 
in Eq.(|9]). The non- vanishing other elements define the numerical accuracy, i.e. the average error 
(6/463 + 6/182 + 6/156)/3 should be around 3%. Since one gets the same 3% numerical error in 
the case of bulk Bl TiN and a negligible one for Bl AIN, we can assume that 3% is numerical 
error threshold for all of our results through the following analysis. One can read from the data in 
Appendix B that some of the SQS supercells result in large non-cubic elements and large deviations 
between the principal cubic elastic constants, which means a breakdown of the cubic symmetry 
relations in Eq.(|9]). Nevertheless, by the previously introduced projection we can extract the closest 
cubic elastic tensors and calculate the distance variations \\E — E'^^'^-\\j\\E\\. These deviations are 
shown in FigjlJ The required 9 elastic constants are summarized in Table III, while the obtained 



projected cubic elastic constants are listed in Table IV Since the C3-(2 x 2 x 2) supercell should 
have cubic symmetry, its \ \E — i?'^^^^'||/||£'|| value defines the numerical threshold for the deviations, 
which is around 4.3%. Thus, only the (2 x 2 x 2), C3-(2 x 2 x 2), Bl-(2 x 2 x 2) and (4 x 4 x 3) 
supercells give cubic symmetry within the most general 21 dimensional vector space related to the 
21 elastic constants. The (4 x 3 x 2)* and (4x4x4) are the candidates to exhibit closely to cubic 
symmetry from elastic point of view. The ordered LIq structure results in the largest deviation 
from cubic symmetry. While the (2x2x2) supercell with relative large SRO parameters fulfills 
the cubic requirement, the larger and perfectly random (4x3x4) supercell does not. In general it 
underlines the importance of applying supercells designed with the inclusion of symmetry in modeling 
anisotropic tensorial properties of alloys. Namely, the SRO parameters or the atomic configuration 
should be optimized in such a way as to support also the point group symmetry. From Fig. [T] with 
including the SRO parameters we conclude that among the tested supercell structures our (4x4x3) 
model should be taken as the closest SQS model to study the elasticity in cubic Tio.5Alo.5N. Thus, 
we conclude that the accurate elastic constants of Tio.5Alo.5N alloy are C\\ = 447 GPa, C12 = 158 
GPa and C44 = 203 GPa within 3% of numerical error. We also see, that with using very ad-hoc or 
inadequate structures, such as the LIO, in predicting elastic constants of Tio.5Alo.5N one faces with 
large 22-50% errors. 

The projection technique allows us to evaluate the supercells in a smaller, 9 dimensional vector 
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space. In the following we consider only the nine elastic constants of Cn, C22, C*33, C12, C13, C23, 
C44, C55 and Cqq. These elastic constants are given in Table [111} The deviations of these constants 
from the projected cubic elastic constants are shown in Fig. [2j In this figure the three columns 
for each supercell give the deviations along the three orthogonal directions, [100], [010] and [001]. 
One can see in the figure, where the horizontal lines show our 3% error threshold, that in the 9 
dimensional space only three supercells, the (2 x 2 x 2), C3-(2 x 2x) and (4x4x4) give cubic 
elastic symmetry. Similarly to Fig. [l] the (4x4x3) supercells performs very well, while the totally 
random (4x3x4) does not. Accordingly, Fig. |2] correlates quite well with Fig. [T| namely we see 
the same set of structures that performing perfectly good or bad. This leads us to the conclusion 
that one can analyze the performance of the supercells in describing cubic elasticity within this 9 
dimensional subspace, too. This means a great reduction in the computational cost, since only 9 
elements have to be calculated to measure the representation of elasticity. By the way, the analysis 
in this 9 dimensional subspace might result in another best approximate superstructure, like in this 
study. Fig. [2] shows clearly, that the (4x4x4) supercell results in a somewhat better representation 
of cubic elastic symmetry in this space. 

However, this small discrepancy between the two previously performed analysis, within the full 21 
and 9 dimensional spaces can be resolved by comparing the derived projected cubic elastic constants 
of the supercells. This comparison is shown in Fig. [3j In Fig. [3] the relative deviations of elastic 
constants are plotted with respect to the values obtained for the (4x4x3) SQS in correspondence 
with the conclusion from Figj2] As one can see, the (4x3x4) and (4x4x4) supercells actually 
result in the same cubic elastic constants within the 3% numerical error. An interesting fact is that 
the values in Fig. [3] should correlate with the corresponding relative differences in Fig. [Tj See, for 
example, the big difference between the cases of (4 x 4 x 4) and Bl-(2 x 2 x 2). Accordingly, Fig. 
[3] concludes the convergency of the cubic elastic constants of Tio.5Alo.5N with respect to differently 
shaped and sized supercell models. Accordingly, the projected cubic elastic constants can be used to 
predict elasticity of cubic alloys. 

Since the elastic anisotropy in TiAlN alloys has a huge impact on the materials mechanical prop- 
erties [7], an accurate prediction of the Zener's elastic anisotropy is of a big importance. Using the 
projected cubic elastic constants one can derive the Zener's elastic anisotropy via 

^=7.^. (10) 

The derived values are listed in Table IV and plotted in Fig. |4} The 3% numerical error accumulates 
in the nominator and results in the approximate 5% difference between the two elastic anisotropy 
values obtained with the best (4x4x4) and (4x4x3) supercells. Thus, the Zener's elastic 
anisotropy in Tio.5Alo.5N should have the value of A=1.40 with around 6% numerical error. Fig|4] 
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shows not only the cubic projected elastic anisotropy values but also their variation along the three 



orthogonal directions, [100], [010] and [001], using the data from Table III, For example, in the [100] 
direction one has y4(x, y, 2) = 2C44/(Cii — C12) while in the [010], yl(?/za;) = 2C66/(C'33 — C13). These 
orientational variations should vanish in case of true cubic symmetry. However, as the figure shows 
one may get a large orientation dependence (~ 55%, see Cl-(2 x 2 x 2)) for a supercell being far 
from fulfilling cubic point symmetry. The sizes of the variations should correlate with the deviations 
shown in Fig. [2j Accordingly, Fig. |4] gives a similar way to analyze the performance of the supercells 
in modeling elasticity of cubic systems. 

In the Reuss averaging method, with assumed uniform stress distribution, the strain ratio 
e[2oo]/e[iii] = -E'[iii]/-£'[ioo]) where E[hki] denotes the directional Young's elastic moduh, can be ap- 
plied to estimate elastic anisotropy experimentally in cubic materials. Using our most accurate 
supercell model of (4 x 4 x 3) the strain ratio £[200] /^[iii] has the value of 1.32 in Tio.5Alo.5N. It is 
also shown in Fig|4| This value deviates from our A = 1.40 value less than the 6% numerical error. 
The elasticity of polycrystalline Tio.5Alo.5N can be discussed in terms of the Reuss and Voigt bulk 
(i?R, By) and shear moduli (Gr, Gy) and also the derived Young's modulus (_Ev/g) and Poisson ratio 
(j^v/g), 
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V 



^ (Cn + G22 + G33) + 2(Ci2 + Ci3 + C23) 

= 9 

(Cii + C22 + C33) - (C12 + Ci3 + C23) + 3(C44 + C55 + Ces) 

15 



(5*11 + 5*22 + >S'33) + 2(5*12 + Sis + 'S'23) 

^ ^ 15 

4(^n + ^22 + 533) - 4(^12 + 5i3 + 523) + 3(^44 + ^55 + See) ' 

9BG W-2G 

h = , v= 11 

W + G 6B + 2G ^ ' 

where Sij denotes the elastic compliances. These polycrystalline averaged quantities obtained for our 
the supercell (4x4x3), that approximates both the randomness and cubic symmetry accurately, 
are summarized in Table IVl The values clearly show the cubic requirement of -Br = By. 



SUMMARY 



In this study we discuss the performance of superstructures, including approximate special quasir- 
andom structure (SQS) supercells in modeling the elasticity of cubic Bl Tio.5Alo.5N alloy. Though 
the SQS approach provides a successful scheme to model and predict the thermodynamics of al- 
loys, the technique is not aiming to represent tensorial materials properties with symmetry. Thus, 
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its straightforward application can not provide an unambiguous description of elasticity in random 
alloys. 

Here, we applied a symmetry based projection technique to accurately predict the cubic elastic 
tensor of Bl Tio.5Alo.5N alloy within the SQS approach. We derived from ab-initio calculations the 
closest cubic elastic tensor of Bl Tio.5Alo.5N by using several supercells. With the help of these derived 
cubic projected elastic constants we presented a detailed analysis and comparison of differently 
shaped and sized supercell models in describing elasticity of a system with cubic symmetry. Thus, 
we accurately determined the cubic elastic constants of cubic Tio.5Alo.5N. The (4x4x3) supercell 
provided us the best model of both, randomness and elasticity, which resulted in Cu = 447 GPa, 
C12 = 158 GPa and C44 — 203 GPa for the cubic elastic constants with 3% of error and A — 1.40 for 
Zener's elastic anisotropy with 6% of error. 

With the help of the obtained elastic tensors, each with 21 constants, our results established the 
fact that supercells with good SRO parameters may include large non-cubic elastic constants and, 
on the contrary, supercells with bad SRO parameters might approximate cubic elastic tensor fairly 
accurately. We showed that using only 9 elements, Cu, C22, ^33, C12, C13, C23, C44, C55 and Cee 
constants from the tensors, one can also adequately evaluate the supercell models and convergency of 
the results. We also showed that, the deviations between the three equivalent Zener-type anisotropy 
factors, oriented along the [100], [010] and [001] directions, confirm the same observation and establish 
a measure of approximate cubic symmetry. 

In summary, in this study we accurately predict cubic elastic constants of Bl Tio.5Alo.5N alloy 
and establish in general the importance of selecting proper SQS supercells with symmetry arguments 
to reliably model elasticity of alloys. Furthermore, we suggest the calculation of nine elastic tensor 
elements - Cu, C22, C33, C12, C13, C23, C44, C55 and Cee, to evaluate and analyse the performance 
of supercells in describing elasticity of alloys. 
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Appendix A: The atomic distributions and coordinates in the supercells 



The atomic distributions and coordinates relative to the supercells lattice parameters are listed 
in Tables Vl|VII For the atomic distributions and coordinates in Cl-(2 x 2 x 2), C3-(2 x 2 x 2) and 
Bl-(2 X 2 X 2), see Ref.[20]. 
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Appendix B: Elastic tensors (in GPa) of cubic TiN, AIN and Tio.5Alo.5N calculated for 
supercells from Table |l] 



Elastic tensor of Bl TiN 



Elastic tensor of Bl AIN 



Elastic tensor of LIq structure 



( 617 


123 


123 


-2 


-6 


-6\ 
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-6 
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-6 


-2 
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Elastic tensor o/ (2 x 2 x 2) SQS 



Elastic tensor o/ (2 x 3 x 2) SQS 



Elastic tensor o/ (4 x 3 x 2) SQS 



^469 


148 


151 


-3 


-5 


-3\ 




488 


148 


-3 
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-3 






469 


-3 


-5 


-3 
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-4 


-4 










208 
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210y 


^429 


173 


164 


2 


4 


6\ 




388 


169 


15 


16 


4 






443 


11 


9 


16 
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9 


8 










203 


9 
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188y 


^436 


161 


160 


12 


11 


25 \ 




453 


160 


4 


15 


1 






428 


13 


3 


8 








188 


12 


9 










186 


9 


\ 










189y 
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Elastic tensor o/ (4 x 3 x 2)* SQS 

U77 144 155 -2 2 9 \ 

445 149 -1 3 -14 

474 3 -6 1 

210 2 

215 1 

^ 199 J 

Elastic tensor of Cl-{2 x 2 x 2) structure 

^385 164 164 4 4 4 \ 

495 136 1 

495 1 

222 3 3 

183 2 

V 183; 

Elastic tensor of C3-{2 x 2 x 2) structure 

/462 156 156 6 6 6 \ 

462 156 6 6 6 

462 6 6 6 

182 6 6 

182 6 
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Elastic tensor of Bl-{2 x 2 x 2) structure from Ref.f2U^. 



Asi 139 147 -1 -4 -2\ 

482 147 -4 -1 -2 

473 -1 -1 -2 

214 -1 -1 

214 -1 

\ 201/ 

Elastic tensor o/ (4 x 3 x 4) SQS 



/431 


148 


153 


-2 


25 


21 \ 




478 


148 


-11 


-10 


-16 






472 


9 


-11 


5 








216 


3 













196 


-1 


v 










194 y 



Elastic tensor 0/ (4 x 4 x 3) SQS 

/456 161.05 152 1 2 4 \ 

425 160 7 4 1 

460 3 1 9 

201 3 5 

211 3 

^ 198y 
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Elastic tensor 0/ (4 x 4 x 4) SQS 



/457 149 156 -2 14 19 \ 



462 156 



11 



3 



16 



444 17 



5 



1 



202 



1 



203 



1 



V 



200 / 
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FIG. 1: (Color online) Calculated euclidian norm deviations 1 1 £^ — E"^*^*^ | 
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FIG. 2: (Color online) Comparison of the calculated elastic tensor elements with the projected principal 
cubic elastic constants in Tio.5Alo.5N. 
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FIG. 3: (Color online) The calculated projected cubic elastic constants of Tio.5Alo.5N relative to the values 
obtained for the (4x4x3) SQS model. 
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FIG. 4: (Color online) Zener's elastic anisotropy values in Tio.5Alo.5N for each structural models considered 
in this study. The horizontal solid line shows the value of A obtained for the (4x4x3) SQS. 
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TABLE L The Warren-Cowley pair short range order parameters (SROs) up to the 7th neighboring shell 
for each SQS supercell considered in this work. 



str.\shell 


number of atoms 


1 


2 


3 


4 


5 


6 


7 


Llo 


8 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 


(2x2x2) 


16 


-0.16 


0.0 


-0.16 


1.0 


-0.16 


0.0 


-0.16 


(2x3x2) 


24 


-0.11 


0.0 


-0.08 


0.33 


-0.06 


-0.08 


0.03 


(4x3x2) 


48 


0.0 


0.0 


0.0 


0.0 


-0.06 


0.0 


0.0 


(4x3x2)*'^ 


48 


-0.14 


0.28 


-0.10 


0.14 


0.0 


-0.17 


-0.01 


Cl-(2x2x2) 


64 


-0.33 


1.0 


-0.33 


1.0 


-0.33 


1.0 


-0.33 


C3-(2x2x2) 


64 


0.0 


-1.0 


0.0 


-1.0 


0.0 


-1.0 


0.0 


Bl-(2 X 2 X 2)^ 


64 


0.0 


0.0 


0.0 


-0.33 


0.0 


0.0 


0.0 


(4x3x4) 


96 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


(4x4x3) 


96 


0.0 


0.0 


-0.01 


-0.01 


-0.01 


0.0 


0.0 


(4x4x4) 


128 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



"I'he * marks a ditferent atomic conhguration in the supercell. 
''The supercell was obtained by von Pezold et al. in Ref. [20] . 



TABLE II: The optimized structural parameters of Tio.5Alo.5N obtained with different supercells from Table 



str. 


a(A) 


6(A) 


c(A) 




Z(a, c) 


l{a,b) 


Llo 


4.17 


4.17 


4.24 


90.00 


90.00 


90.00 


(2x2x2) 


4.18 


4.18 


4.22 


59.72 


59.72 


59.73 


(2x3x2) 


4.20 


4.19 


4.17 


60.17 


60.12 


59.65 


(4x3x2) 


4.18 


4.18 


4.20 


59.67 


59.79 


59.99 


(4x3x2)*'^ 


4.19 


4.19 


4.17 


60.14 


60.17 


59.66 


Cl-(2x2x2) 


4.25 


4.16 


4.16 


90.00 


90.00 


90.00 


C3-(2x2x2) 


4.18 


4.18 


4.18 


90.00 


90.00 


90.00 


Bl-(2x2x2)'' 


4.18 


4.18 


4.18 


89.78 


89.78 


90.00 


(4x3x4) 


4.20 


4.15 


4.18 


60.26 


59.91 


60.38 


(4x4x3) 


4.19 


4.18 


4.16 


60.15 


60.13 


59.94 


(4x4x4) 


4.18 


4.18 


4.19 


59.93 


60.08 


60.00 



"The * marks a ditterent atomic conhguration in the supercell. 
''The supercell was obtained by von Pezold et al. in Ref. [50] . 
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TABLE IIL The calculated GGA elastic tensor elements of Tio.5Alo.5N from Appendix B, whichs are involved 



in the cubic project: 


ion, see 


Eqs.Q and 














str.\const. 


Cn 


C22 


C33 


C12 


Cl3 


C23 


C44 


C55 


Cee 


TiN 


617 


618 


618 


123 


123 


123 


178 


178 


178 


AIN 


402 


402 


402 


157 


157 


157 


300 


300 


300 


Llo 


409 


409 


332 


183 


197 


197 


100 


100 


120 


(2x2x2) 


469 


488 


469 


148 


151 


148 


210 


208 


210 


(2x3x2) 


429 


388 


443 


173 


164 


169 


187 


203 


188 


(4x3x2) 


436 


453 


428 


161 


160 


160 


188 


186 


189 


(4x3x2)*'' 


477 


445 


474 


144 


155 


149 


210 


215 


199 


Cl-(2x2x2) 


385 


495 


495 


164 


164 


136 


222 


183 


183 


C3-(2x2x2) 


462 


462 


462 


156 


156 


156 


182 


182 


182 


Bl-(2x2x2)'' 


481 


482 


473 


139 


147 


147 


214 


214 


218 


(4x3x4) 


431 


478 


472 


148 


153 


148 


216 


196 


194 


(4x4x3) 


456 


425 


460 


161 


152 


160 


201 


211 


198 


(4x4x4) 


457 


462 


444 


149 


156 


156 


202 


203 


200 



"The * marks a ditferent atomic configuration in the supercell. 
''The supercell was obtained by von Pezold et al. in Ref. [20] . 



TABLE IV: The projected principal cubic elastic constants and the derived Zener's elastic anisotropy A = 



2C44/ (Cn — C12) 


of Tio.5Alo.5N obtained with the different structural 


models. 




str.\const. 


Cn 


C12 


C44 


A 


Llo 


384 


193 


107 


1.12 


(2x2x2) 


475 


149 


209 


1.28 


(2x3x2) 


420 


169 


193 


1.53 


(4x3x2) 


439 


160 


188 


1.35 


(4x3x2)*'' 


465 


149 


208 


1.32 


Cl-(2x2x2) 


459 


155 


196 


1.29 


C3-(2x2x2) 


462 


156 


182 


1.19 


Bl-(2x2x2)^ 


479 


144 


215 


1.29 


(4x3x4) 


460 


150 


202 


1.30 


(4x4x3) 


447 


158 


203 


1.40 


(4x4x4) 


454 


154 


202 


1.34 
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"The * marks a different atomic configuration in the supercell. 
''The supercell was obtained by von Pezold et al. in Ref. [20] . 



TABLE V: The polycrystalline bulk (B), shear (G), Young (E) moduli in unit of GPa and the Poisson ratio 



of Ti, 


0.5AI0.5N obtained with the (4x4x3) supercell. 










Bv 


Br 


Gv 


Gr 


Ey 


Er 






254 


254 


180 


174 


437 


425 


0.21 


0.22 


TABLE VL The internal atomic structures of the supercells in relative coordinates. 




I Z A Z A Z ) 




l'9 V V 9^ 

I Z A A Z I 




/'zl V V 9^1 
1 ^ A A Z ) 
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Ti 
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n n n 
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Al 





Al 


1 /9 n n 


Al 
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Al 
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i/ 4 U U 


Ti 
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Ti 
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n 1 /9 n 


Ti 
X 1 


n 1 /*? n 
u 1/ u 


Al 
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Al 
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Ti 


1 /9 1 /9 n 
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Al 


n 9 /Q n 

U Z/ u 


Ti 
X 1 


q //I n n 
0/ 4 U U 


Ti 
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Al 
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Ti 
J. 1 


1 /9 1 /Q n 
i/ Z i / u 
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Ti 
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Ti 


1 /9 n 1 /9 
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Ti 
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Al 
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Al 


n n 1 /9 

U U i/ z 


Ti 
X 1 


1 //I 1 /"? n 
i/ 41: i/ u 


Al 


1/4 1/3 


Al 


1/91/91/9 
i/Z 1/Z 1/Z 


Al 


1 /9 n 1 /9 

1/ Z U i/Z 
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1 //I 9 /*? n 

i/ 4 Z/ u 


Ti 
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Ti 
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n 1 /Q 1 /9 

U i/ i/Z 
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1 /9 1 /Q n 
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Ti 
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A 1 


q //I 1 /q n 

0/4 i/O U 
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Ti 3/4 2/3 1/2 Al 3/4 2/3 1/2 



TABLE VII: The internal atomic structures of the supercells in relative coordinates (continuation). 



(4x3x4) (4x4x3) (4x4x4) 



Ti 





Al 





Ti 





Al 


1/4 


Al 


1/4 


Al 


1/4 


Ti 


1/2 


Al 


3/4 


Al 


1/2 


Ti 


3/4 


Al 
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Ti 
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Ti 
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Al 
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Ti 
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Ti 
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Al 
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Ti 
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Al 


1/4 1/3 


Al 
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Al 
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Al 
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Al 
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